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The electric dipole {El) strength functions have been systematically calculated 
based on the time-dependent density functional theory (TDDFT), using the 
finite amplitude method and the real-time approach to the TDDFT with pair- 
ing correlations. The low-energy El strengths in neutron-rich isotopes show 
peculiar behaviors, such as sudden enhancement and reduction, as functions 
of the neutron numbers. They seem to be due to the interplay between the 
neutron shell effect and the deformation effect. 
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1. Introduction 

The electric dipole {El) response is a fundamental mode of excitation in 
nuclei and a useful tool to probe the isovector property of nuclei. The giant 
dipole resonance may provide information on the symmetry energy near 
the saturation density po- In contrast, the low-energy El modes, which are 
often referred to as pygmy dipole resonances, is sensitive to the nuclear 
structure and may carry information on the symmetry energy £^sym(p) at 
densities away from pQ. Thus, the low-energy dipole modes in neutron-rich 
nuclei are of significant interest at present. 

The low-energy dipole strengths have been experimentally observed in 
several neutron-rich isotopes; O isotopes, ^'^ ^^Ne,'^ and Sn isotopes. The 
low-energy El strength observed in stable isotopes carries only less than 
< 1 % of the Thomas-Reiche-Kuhn (TRK) sum-rule value. In contrast, it 
may amount up to about 5 % in neutron-rich nuclei. Therefore, we may 
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expect a qualitative difference in properties of the low-energy dipole modes 
of excitation between in stable and in neutron-rich nuclei. 

The time-dependent density-functional theory (TDDFT) provides us 
with a practical and reasonable description of nuclear strength functions.^ 
Its rigorous theoretical foundation is given by the one-to-one correspon- 
dence between the time-dependent external potential v{t) and the time- 
dependent one-body density p{t).^^ In the Kohn-Sham scheme, it gives the 
time-dependent Kohn-Sham (TDKS) equations that are practically identi- 
cal to the one known as the time-dependent Hartree-Fock (TDHF) equa- 
tions in nuclear physics. An extension for the study of superconducting 
systems has been also carried out by including the time-dependent pair po- 
tential A(t) and the time-dependent pair density This eventually 
leads to the "time-dependent Bogoliubov-de-Gennes-Kohn-Sham (TDB- 
dGKS) scheme" , which is, in nuclear physics, much more familiar with the 
name of "time-dependent Hartrce-Fock-Bogoliubov (TDHFB) equation". 
The TDBdGKS equations determine the time evolution of quasiparticle or- 
bitals whose number is 2M where M indicates the dimension of the single- 
particle model space. Since M is significantly larger than the particle num- 
ber N, it requires a huge computational task and becomes a challenging 
subject in computational nuclear physics. ^'^^ 

In this paper, alternative approaches based on the TDDFT, which are 
computationally more feasible than the full solution of the TDBdGKS equa- 
tions, are utilized to study the properties of the low-energy El strength. 
These methods will be briefly explained in Sec. 2. The numerical applica- 
tions are given in Sec. 3. 

2. Theoretical tools 

We use the canonical-basis formulation of the TDBdGKS method and the 
finite amplitude method (FAM) to obtain low-energy El strength func- 
tions in neutron-rich nuclei. We briefly recapitulate these methods in the 
foUowings. Details of the method can be found in references given below. 

2.1. Canonical-basis TDBdGKS equations 

The real-time calculations of the TDKS equations have been carried out 
in real space. ^^'^^'^^ However, the computational cost is significantly in- 
creased by inclusion of the pairing correlation (TDKS — > TDBdGKS), which 
makes practical applications very difficult. This numerical cost can be re- 
duced, by several orders of magnitude, introducing an approximation for 
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the time-dependent pair potential, similar to the BCS approximation in 
static cases. This may lead to the following set of equations: 



i§i\Mt)) 


= (/i(i)-%(i))|<^fe(i)>, 


(1) 






(2) 






(3) 




= {rjkit) + rjkit)) Kk{t) + Afe(t) {2pk{t) - 1) . 


(4) 



Here, Eqs. (1) and (2) describes the time evolution of a pair of canonical 
states, \(l)k{t)) and (t)). Their occupation and pair probabilities are given 
by Eqs. (3) and (4). The Hamiltonian h{t) is a functional of density h[p{t)], 
and r/fc {t) and r/j {t) are arbitrary real parameters. In numerical calculations 
in Sec. 3, we adopt r]k{t) = {Mt)\hit)\Mt)) Vkit) = ('/'fe Wl/^WI^fe W)- The 
computational task for solution of these equations is roughly similar to that 
of TDKS equation. 



2.2. Finite amplitude method for linear response 

The finite amplitude method (FAM)^^ is another feasible approach to the 
linear response calculation. The formulation has been extended to super- 
fluid systems as well.^® The method allows us to easily construct a computer 
code for the linear response calculation based on the TDDFT. The FAM 
has been applied to the coordinate-space representation, -'^^"^-'^ and to the 
harmonic-oscillator-basis representation.^^ The essential idea is that the 
complicated residual induced part of the Hamiltonian 5h can be calculated 
in terms of the finite difference associated with the non-Hermitian density 
prj. Provided that \(t>k) are the canonical single-particle states at the ground 
state {ho\(j)k) = ek\(l)k) and po = I]i:hoio l^l^i) {<l>i\) ^ 

5h = ^ {h[pr,] - ho) , (5) 
Pri = Po + rjdp = ^ l^iXV'il, (6) 



z:hole 



where 



\^i) = \4>i)+ri ^rnMm), (7) 

m:particle 

{^i\ = {(j)i\+r] "^mi{<i>m\- (8) 

m:particle 
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In the FAM, the induced field Sh for a given (X, Y) can be estimated using 
Eq. (5). This only requires us to calculate the single- particle (Kohn-Sham) 
Hamiltonian h[p] with different bra's and ket's in Eq. (6). Then, we resort 
either to iterative algorithms^'"'^^'^^ or to the matrix diagonalization,^-^ in 
order to obtain solutions of the linear- response equation {X,Y). 

3. Low-energy El strength 

The El strength functions in even-even isotopes are systematically calcu- 
lated. For relatively light nuclei, the FAM was used neglecting the pairing 
correlations. ^'^ For heavier nuclei, in which the pairing is expected to play 
an important role, we use the real-time method based on the canonical-basis 
TDBdGKS method with a time-dependent perturbation of the external El 
field.15'16 

3.1. Shell effects 

The photoabsorption cross sections for Ca isotopes (N = 20—40), estimated 
in the dipole approximation, are shown in Fig. 1. For Ca isotopes with 
A'' = 20 — 28, we see only neghgible El strengths below E = 10 MeV. In 
contrast, those with N = 30 — 40 shows sizable El strengths below 10 McV. 
This indicates that a sudden appearance of the low-energy El modes takes 
place beyond the magic number A^ = 28 in Ca isotopes. We have confirmed 
the same behavior in neighboring isotopes. 

Similar jumps of the low-energy El strengths can be observed at A^ = 
14 — > 16 and A' = 50 — )■ 52.^° These indicate a strong shell effect on the 
appearance of the low-energy El strength. These numbers are related to 
the occupation of single-particle orbitals with low orbital angular momenta, 
such as Si/2 (A^ = 14 16), P3/2 {N = 28 ^ 30), and ^5/3 (A^ = 50 -> 
52). When these \ow-£ orbitals are weakly bound having spatially extended 
characters, we may expect the threshold effect which may enhance the El 
strengths near the neutron emission threshold energy. Thus, we suppose 
that the low-energy El strengths predominantly possess a single-particle 
nature. 

3.2. Effects of deformation and separation energy 

Let us move toward heavier isotopes. We use the canonical-basis TDBdGKS 

method to calculate the El strength distribution.^'^® We have found that 
the effect of pairing correlation is not so significant, in general. However, 
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Fig. 1. Photoabsorption cross sections in Ca isotopes, calculated using the FAM with 
the SkM* functional. The origin of the vertical axis is shifted upward by 30 mb as the 
neutron number increases by two. 



for selected nuclei, the pairing correlation affects the shape of the ground 
state, which modifies the El properties accordingly. 

We calculate the integrated energy-weighted El strength up to _E = Ec, 
defined as 

In Fig. 2, the ratio of toi(IO MeV) to toi(oo) is shown for Sr isotopes. The 
ratio is less than 1 % for isotopes with N < 50. Then, the ratio jumps 
up beyond N = 50, which is consistent with the argument given above. 
However, there are sudden drops of the low-energy El ratio at = 58 — > 60 
and N — 72 —> 74. This seems to be due to changes of the ground-state 
deformation. 
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Fig. 2. mi {i?)/mi (oo) with E = 10 MeV for even-even Sr isotopes. The calculation was 
performed with the SkM* functional using the canonical-basis TDBdGKS equations. See 
text for details. 



The calculated ground states in even-even Sr isotopes with TV = 40 — 50 
are all spherical {jS — 0). The two-neutron separation energies, which are 
equal to twice of the chemical potential, gradually decrease as the neu- 
tron number increases. Then, the calculation predicts that the onset of the 
ground-state deformation takes place at TV = 58 -> 60, from spherical to 
prolate shapes (/3 = — ^ 0.37). This shape transition result in the increase 
of the two-neutron separation energy and the decrease of the low-energy El 
strengths. The deformation stays rather constant for N — 60 — 72, with pro- 
late shapes of /3 = 0.34 — 0.38 and with decreasing separation energies as in- 
creasing the neutron number. At = 72 — > 74, again, the shape transition 
takes place, from prolate to oblate shapes, (/3,7) = (0.34,0) (0.14,60°). 
This shape change accompanies the increase of the separation energy and 
decrease of the low-energy El strength. 

The low-energy strength can be decomposed into those associated with 
X, y, and z directions, as in the last equation in Eq. (9). This decomposition 
is also shown in Fig. 2. The z axis is chosen as the symmetry axis for axially 
deformed nuclei. For prolate Sr isotopes with iV = 60 — 72, the strength 
associated with the z (K — 0) component is dominant. The K ^ domi- 
nance was also reported for neutron-rich Sn isotopes using the relativistic 
quasiparticle random phase approximation.^"^ This was interpreted by the 
conjecture that the neutron skin is thicker in the z direction than the x and 
y directions. We calculate the neutron-skin thickness in the x direction as 
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^Xnp = sj {x'^)n — \/ (a;^)p, and those for y and z directions in exactly the 
same way. It turns out that Ar„p is even larger with respect to the x (y) 
direction than the z direction, for prolate nuclei. Therefore, the observed 
K = Q dominance in the low-energy El strength cannot be attributed to 
the different neutron skin thickness. This suggests that these low-energy 
strengths are not associated with the skin modes. 

4. Summary 

The El strength functions have been systematically calculated with the 

finite amplitude method and the real-time method, based on the time- 
dependent density functional theory. The low-energy El strength distri- 
butions in stable and neutron-rich isotopes were estimated from these cal- 
culations. We have found a strong neutron shell effect and have identified 
magic numbers for the appearance of low-energy El modes. The deforma- 
tion and separation energies also play an important role in the low-energy 
El strength distribution. 
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